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Weak radiative decays of the B mesons belong to the most important flavor changing processes 
that provide constraints on physics at the TeV scale. In the derivation of such constraints, ac¬ 
curate standard model predictions for the inclusive branching ratios play a crucial role. In the 
current Letter we present an update of these predictions, incorporating all our results for the C>(o! 2 ) 
and lower-order perturbative corrections that have been calculated after 2006. New estimates of 
nonperturbative effects are taken into account, too. For the CP- and isospin-averaged branching 
ratios, we find B s 7 = (3.36 ± 0.23) x 10“ 4 and Bd 7 = {lN^ho' 22 ) x 10 —5 , for E y > 1.6 GeV. Both 
results remain in agreement with the current experimental averages. Normalizing their sum to the 
inclusive semileptonic branching ratio, we obtain R-, = (B S7 + Bd-y) IB c i v = (3.31 ± 0.22) x 1CP 3 . 

A new bound from B s -y on the charged Higgs boson mass in the two-Higgs-doublet model II reads 
M h ± > 480 GeV at 95%C.L. 

PACS numbers: 13.20.He, 12.38.Bx, 12.60Fr 


I. INTRODUCTION 

The inclusive decays B —» X s j and B — > X^pf are con¬ 
sidered among the most interesting flavor changing neu¬ 
tral current processes. They contribute in a significant 
manner to current bounds on masses and interactions of 
possible additional Higgs bosons and/or supersymmetric 
particles. The evaluation of such bounds depends in a 
crucial manner on both the central values and uncertain¬ 
ties of the branching ratio predictions within the stan¬ 
dard model (SM). Updating the SM predictions is the 
main purpose of the present Letter. 

Measurements of the CP- and isospin-averaged B —>• 
X s "f branching ratio by CLEO [l|, Belle j| |f|, and 
BABAR j(i lead to the combined result [dj] 

S® 7 P = (3.43 ± 0.21 ± 0.07) x 10 -4 , (1) 

for the photon energy Ej > Eo = 1.6 GeV in the de¬ 
caying meson rest frame. The combination involves an 
extrapolation from measurements performed at Eq £ 
[1.7, 2.0] GeV. Applying the same extrapolation method 
to the available B -a X^P) measurement [l0|, one finds 

( 2 ) 


at E 0 = 1.6 GeV E2- More precise determinations of 
U® 7 P for q = s,d are expected from Belle II fl2| . 

Theoretical calculations of B ql have a chance to match 
the experimental precision only in a certain range of Eq 
where the nonperturbative contribution 5T nonp in the re¬ 
lation 

T(H -I X ql ) = T(b -1 X p 7 ) + sr nonp (3) 

remains under control. Here, T(fe —»• X p 7 ) denotes the 
perturbatively calculable rate of the radiative 6-quark de¬ 
cay involving only charmless partons in the final state. 
Their overall strangeness vanishes for X 3 ! and equals —1 
for Xf. The analysis of Ref. [l3| implies that unknown 
contributions to <5r nonp are potentially larger than the so- 
far determined ones, and induce around ±5% uncertainty 
in B S1 at Eq = 1.6 GeV. Nonperturbative uncertainties in 
Bdj receive additional sizeable contributions [l|J due to 
collinear photon emission in the b -P duup process whose 
Cabibbo-Kobayashi-Maskawa (CKM) factor is only a few 
times smaller than the one in the leading term. 

Apart from possible future progress in analyzing non¬ 
perturbative effects, one needs to determine T(6 —> A p 7 ) 
to a few percent accuracy. It requires evaluating next- 
to-next-to-leading order (NNLO) QCD corrections that 
involve Feynman diagrams up to four loops. The first 


p = (1.41 ± 0.57) x 10 -5 
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SM estimate of the B —> X s 7 branching ratio at this level 
was presented in Ref. [15] where all the corrections calcu¬ 
lated up to 2006 were taken into account. A part of the 
0(a^) contribution was obtained via interpolation [l6j] 
in the charm quark mass between the large-m c asymp¬ 
totic expression 0 and the m c = 0 boundary condition 
that was estimated using the Brodsky-Lepage-Mackenzie 
(BLM) approximation [18 |. 

In the present Letter, we provide an updated predic¬ 
tion for £> S7 , including all the contributions and estimates 
worked out after the completion of Ref. 0 [. They are 
listed in Sec. El where the necessary definitions are in¬ 
troduced. The interpolation in m c is still being applied. 
However, the m c = 0 boundary condition is no longer a 
BLM-based estimate but rather comes from an explicit 
calculation 0- 

The current analysis supersedes our previous one in 
Ref. 0, which was published in 2006 and has not been 
updated since then. It has been widely considered as 
a standard reference until now. The time for our up¬ 
date comes only at present because the most recent and 
technically challenging four-loop calculation of Ref. [l9|] 
constitutes a breakthrough in the analysis. It has an im¬ 
portant effect on the central value of B sl . 

The Letter is organized as follows. After discussing 
B S1 in Sec. El our NNLO analysis is extended to 
Bd-y in Sec. Mil Next, in Sec. [[Vi we consider i? 7 = 
(B S1 ± Bd-y) IBctv which may sometimes be more conve¬ 
nient than B SJ for deriving constraints on new physics. 
Sec. [V] is devoted to presenting a generic expression for 
beyond-SM contributions, as well as an updated bound 
for the charged Higgs boson mass in the two-Higgs- 
doublet model II (THDM II). We conclude in Sec. IVII 


II. B S7 IN THE SM 


Radiative R-meson decays are most conveniently de¬ 
scribed in the framework of an effective theory that arises 
after decoupling of the W boson and heavier particles. 
Flavor-changing weak interactions that are relevant for 
T(6 —)• X^y) with q = s,d are given by 


£eff ~ v; q v tb 


CJiQi + Ci{Qi - Qi) 


• ( 4 ) 


Explicit expressions for the current-current (Q i^), four- 
quark penguin (< 33 ,..., 6 )> photonic dipole {Q 7 ), and 
gluonic dipole (Qs) operators can be found, e.g., in 
Eq. (2.5) of Ref. [lj|. The CKM element ratio n q = 
(VuqVub) / iYtqVtb) is small for q = s, and it affects B S1 by 
less than 0.3%. Barring this effect and the higlrer-order 
electroweak ones, T(6 — > Xfy) in the SM is given by a 
quadratic polynomial in the real Wilson coefficients Ci 


r {b X*7) 


E 

*1.3=1 


CiCj Gij. 


( 5 ) 


A series of contributions to the above expression from 
our calculations in Refs. iTiil - [ 28 l makes the current anal¬ 
ysis significantly improved with respect to the one in 
Ref. [15(. In particular, the NNLO Wilson coefficient 
calculation becomes complete after including the four- 
loop anomalous dimensions that describe Qi ,...,6 -± Qs 
mixing under renormalization 01 ]. Effects of the charm 
and bottom quark masses in loo ps on the gluon lines in 
G 77 0 , G 7 S 0 and G(i 2)7 0 ], as well as a com¬ 
plete calculation of G 73 0 |, are now available. Three- 
and four-body final-state contributions to G$s 0 : 0 
and G(i : 2)8 [26| are included in the BLM approximation. 
Four-body final-state contributions involving the penguin 
and Q™ 2 operators are taken into account at the leading 
order (LO) 0 ] and next-to-leading order (NLO) 01 . 
Last but not least, the complete NNLO calculation [19| 
of G 17 and G 27 at m c = 0 is used as a boundary for 
interpolating their unknown parts in m c . 

Following the algorithm described in detail in Ref. 19], 
taking into account new nonperturbative effects fl3l . 123 . 
Ho], as well as the previously omitted parts of the NNLO 
BLM corrections [0, we arrive at the following SM pre¬ 
diction 

= (3.36 ± 0.23) x 10 -4 for E 0 = 1.6 GeV. ( 6 ) 

Individual contributions to the total uncertainty are of 
nonperturbative (±5%), higher-order (±3%), interpola¬ 
tion (±3%) and parametric (±2%) origin. They are com¬ 
bined in quadrature. The parametric one gets reduced 
with respect to Ref. 0 ] , which becomes possible thanks 
to the new semileptonic fits of Ref. 0]. Our input pa¬ 
rameters, their uncertainties and the corresponding cor¬ 
relation matrix can be found in Appendix D of Ref. ITol . 
Since we normalize to the semileptonic branching ratio 
Bct v , our result shows little sensitivity to the 6 -quark 
mass and the CKM angles. The main parametric uncer¬ 
tainty (±1.5%) originates from B c t Vl while the next one 
(±0.75%) comes from a s (Mz). 

As far as the interpolation uncertainty is concerned, 
one might have hoped for its reduction with respect to 
Ref. [l5| after the explicit evaluation of the m c = 0 
boundary [Till ] . Unfortunately, the interpolated parts of 
the 0(a 8 ) contributions to G(i i 2)7 turn out to be sizeable. 
Their effect on B grows from 0 to around 5% when m c 
changes from 0 up to the measured value (see Fig. 4 of 
Ref. [l9l|). In such a situation, we prefer to stay conser¬ 
vative, and retain our interpolation uncertainty estimate 
at the ±3% level. 

For the higher-order uncertainty estimation, it is use¬ 
ful to study how depends on three renormaliza¬ 

tion scales: the matching scale fj ,0 ~ to* at which the 
heavy particles (t, W, Z, H°) are decoupled, the low- 
energy scale ph ~ mb /2 at which the Wilson coefficient 
renormalization group evolution is terminated, and the 
scale // c at which the charm quark mass is renormal¬ 
ized. We vary them in the ranges /j,q € [80,320] GeV 
and G [1.25, 5] GeV, setting the central values to 

Ho = 160 GeV and Hb = He = 2 GeV. The observed scale 





3 


dependence (see Fig. 6 of Ref. [H|) turns out to be quite 
similar to the one in Fig. 2 of Ref. [ll|. Therefore, we 
leave the higher-order uncertainty estimate at the ±3% 
level, i.e., unchanged with respect to Ref. I 1 . 

The nonperturbative uncertainty estimate of ±5% is 
adopted from Ref. 0 without any modification. It 
turns out to be identical to our earlier rough estimate 
in Ref. [l5[. Some comments on possible future suppres¬ 
sion of this uncertainty are given in Sec. IVII 

The central value in Eq. © is considerably higher than 
3.15 x 10 ~ 4 in Ref. 0 ], although the difference between 
the two values does not exceed the previously estimated 
uncertainty. A detailed description of various contribu¬ 
tions to this difference is given in Sec. 4 of Ref. 0 , as 
well as in Table 2 there. 


III. Bd-v IN THE SM 

Extending our NNLO calculation to the Bdy case be¬ 
gins with inserting the proper CKM factors in Eq. 0 . 
Contrary to k s , the ratio Kd is not numerically small. 
Using the CKM fits of Ref. [33j, one finds 

= (0.007±g;81f) +*(—0.4041°;^). (7) 

The small real part implies that the effects of Hd on the 
CP-averaged Bd-y are dominated by those proportional to 
|«d| 2 . In such terms, perturbative two- and three-body fi¬ 
nal state contributions arise only at the NNLO and NLO, 
respectively. They vanish in the m c = m u limit, which 
effectively makes them suppressed by m 2 /m 2 <0.1. In 
consequence, the main re^-effect comes from b —> duu'y at 
the LO, where phase-space suppression is partially com¬ 
pensated by the collinear logarithms. 

In the first (rough) approximation, one evaluates the 
tree-level b —>■ duu 7 diagrams retaining a common light- 
quark mass m q inside the collinear logarithms [26[ . and 
varying rrib/rn q between 10 ~ and 50 ~ ms/m,r 

to estimate the uncertainty. The considered effect varies 
then from 2% to 11% of Bd-y- A more involved analy¬ 
sis with the help of fragmentation functions gives a very 
similar range [141 ] . Including this contribution in our eval¬ 
uation of the entire Bd-y from Eq. 0, we find 

Bff = (1.73+°;^) x 10 ’ 5 for E 0 = 1.6 GeV, ( 8 ) 

where the central value corresponds to rrib/m q = 50. Our 
result is about 12% larger than the one given in Ref. [Tlj 
where the b —» duuj contributions were neglected. The 
uncertainty estimate in Eq. 0 improves with respect to 
Ref. 0 thanks to including the NNLO QCD corrections 
and using the updated CKM fit 0. Interestingly, the 
parametric uncertainty due to the CKM input amounts 
to ±2.5% only. 

The collinear logarithm problem might seem artificial 
because isolated photons are required in the experimental 
signal sample. Unfortunately, requiring photon isolation 
on the perturbative side would necessitate introducing 


an infrared cutoff on the gluon energies, e.g., in the NLO 
corrections to the dominant G 77 term. Without a dedi¬ 
cated analysis (which is beyond the scope of the present 
Letter), it is hard to verify whether such an approach 
would enhance or suppress the uncertainty in Bd-y- 

Another question concerning the \Kd\ -terms is 
whether the off-shell light vector meson conversion to 
photons can be assumed to be included in our overall 
±5% nonperturbative uncertainty. Much smaller effects 
found in the vector-meson-dominance analysis of Ref. [343 ] 
imply that it is likely to be the case. 

IV. THE RATIO Ry 

In the fully inclusive measurements of radiative 
R-meson decays 0 HH 3 , the final hadronic state 
strangeness is not verified. The actually measured 
quantity is B S y ± Bd-y- Next, the result is divided 
by (l + \(V t * d V t b)/(V t *V t b)\ 2 ) to obtain B S1 . To avoid 
such a complication, we provide here our SM predic¬ 
tion for B S y ± Bd-y with all the correlated uncertainties 
properly taken into account. Moreover, we normalize it 
to the CP- and isospin-averaged inclusive semileptonic 
branching ratio B c t, v . In the B S1 case, such a normal¬ 
ization reduces the parametric uncertainty from ± 2 . 0 % 
to {±1.2,—1.4}%. It may also be useful on the experi¬ 
mental side because the inclusive semileptonic events can 
serve for determining the R-meson yield. Proceeding as 
in the previous sections, we obtain for R 7 = 1.6 GeV 

R® m = ( 6 ™ ± B™) /B clv = (3.31 ± 0 . 22 ) x 10“ 3 . (9) 

The relative uncertainties are identical to those in B S1 
(as given below Eq. ©), except for the parametric one 
which amounts to {±1.2, —1.7}% including the effect of 
m-b/mq. The gain in the overall theory uncertainty is 
hardly noticeable, but this may change with the future 
progress in determining the perturbative and nonpertur¬ 
bative corrections. 


V. BEYOND-SM EFFECTS 

In most of the new-physics scenarios considered in the 
literature, beyond-SM effects on B S1 are driven by new 
additive contributions to the Wilson coefficients of the 
dipole operators at the matching scale /io. Denoting such 
contributions by ACV ,8 and setting /i 0 to 160 GeV, we 
find 

B S y x 10 4 = (3.36 ± 0.23)- 8.22 ACV- 1.99 AC's, 

R 7 x 10 3 = (3.31 ± 0.22)- 8.05ACV-1.94 AC's-(10) 

The above expressions are linearized; i.e., it is assumed 
that the quadratic terms in AC 73 are negligible when 
they enter with 0( 1 ) coefficients into the above equa¬ 
tions. If they are not, a detailed analysis of QCD correc¬ 
tions in the considered beyond-SM scenario is necessary. 
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Such an analysis is available in the THDM II (35j for 
which the NLO f3(ii - [38[ and NNLO |39j corrections to 
AC 7i8 are known. They are always negative and remain 
practically independent of the vacuum expectation value 
ratio tan/3 when tan/3 > 2. Sending tan/3 to infinity in 
the expressions for ACV^, we find the following updated 
bounds from B S1 on the charged Higgs boson mass in this 
model 

M h ± > 480 GeV at 95% C.L., 

M h ± > 358 GeV at 99% C.L. (11) 

For tan /3 < 2 the bounds become considerably stronger, 
but at the same time other observables provide competi¬ 
tive limits poj . In the supersymmetric case, in which the 
charged scalar and the neutral pseudoscalar tend to be al¬ 
most degenerate, the current direct search bounds f44l[42j 
exceed 500 GeV for tan/3 > 20. 

VI. SUMMARY 

We presented an updated prediction for B S1 in the SM 
taking into account all the perturbative and nonpertur- 
bative effects worked out after the 2006 publication Hi 
of the first NNLO estimate for this quantity. Our current 
analysis supersedes the one of Ref. [lg. 

Some of the corrections are still interpolated in 

to c , but the m c = 0 boundary condition now comes from 
an explicit calculation. Despite this improvement, the 
interpolation uncertainty cannot be reduced because the 
interpolated correction is sizeable. Future progress re¬ 
quires extending the calculation of G( 1 , 2)7 to arbitrary 
m c , which is considered a difficult but manageable task. 
It would amount to evaluating the same propagator di¬ 
agrams with unitarity cuts as in Ref. [T§|, but for arbi¬ 
trary m c rather than just for m c = 0. Several hundreds 
of four-loop two-scale master integrals would need to be 
calculated. For this purpose, one could numerically solve 
differential equations in the variable The 

necessary boundary conditions at z 1 could be found 
from asymptotic expansions in this limit. Determining 
such boundary conditions involves only three-loop single¬ 
scale propagator integrals. They are likely much simpler 
than the four-loop single-scale ones in Ref. [Til . 

In parallel, one should investigate whether nonper- 
turbative uncertainties can be suppressed by combin¬ 
ing lattice inputs with measurements of observables like 
the CP- or isospin asymmetries in B —>■ X q "f. In the 
analysis of Ref. (li|, nonperturbative effects have been 
parametrized in terms of the so-called subleading shape 
functions, i.e., matrix elements of nonlocal operators be¬ 
tween the 33-meson states at rest. Determining such 


functions directly seems to remain beyond the current 
lattice capabilities. However, constraints on them can 
be derived from matrix elements of local operators, the 
same ones that matter for the extraction of |V/&| from 
Betv [32| . The higher-dimensional operator matrix ele¬ 
ments are practically unconstrained by the data. Any 
lattice estimates of them could help to suppress the non¬ 
perturbative uncertainties in both \V c b\ and B SJ . 

The main outcome of our current analysis is an up¬ 
wards shift by around 6.4% in the central value of B^. 
It originates mainly from fixing the m c = 0 boundary 
(+3%) and including the complete NNLO BLM correc¬ 
tions to the three- and four-body final state channels 
(+2%). Both effects are within the previously [l5| es¬ 
timated interpolation (±3%) and higher-order (±3%) 
uncertainties, respectively. Nevertheless, the obtained 
0(lcr) increase of the central value is an important one, 
especially in the context of constraining beyond-SM theo¬ 
ries. The new four-loop calculation of the m c = 0 bound¬ 
ary in Ref. 0 improves an essential point in the analy¬ 
sis, and brings the estimated NNLO effects under much 
better control. 

Since is now closer to H“ p (but still B s ™ < B e * p ), 
the bound on M H ± in the THDM II becomes significantly 
stronger. The 95%C.L. one grows by 120 GeV with re¬ 
spect to its previous evaluation in Ref. ^39] (cf. “note 
added” there). For moderate values of tan/3, no other 
available measurement constrains M H ± in a more effi¬ 
cient manner. 

We supplemented our analysis with new NNLO predic¬ 
tions for Bd~f and for the ratio i? 7 = (£> ST + Bdj) IB c i v 
where correlated uncertainties are treated in a consis¬ 
tent manner. The ratio i? 7 may serve in the future as a 
more convenient observable for testing beyond-SM theo¬ 
ries with minimal flavor violation. 
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